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We report the measurement of the branching fractions of the rare decays B° — > Ds*^ + tt~, 
B° -> D ( s * )+ p-, and B {> -> £)£" )_ K (t)+ in a sample of 381 x 10 6 Y(4S) decays into BB pairs collected 
with the BABAR detector at the PEP-II asymmetric-energy e + e~ storage ring. We present evidence 
for the decay B° — > K* + and the vector-vector decays B° — ► D* 3 + p~ and B° — > D* s ~ K* + , as 
well as the first measurement of the vector meson polarization in these decays. We also determine 
the ratios of the CKM-suppressed to CKM-favored amplitudes r(D^n) and r(D l -*'p) in decays 
B° — > Z)'*' ± 7r =F and B° — > D^^p*, and comment on the prospects for measuring the CP observ- 
able sin(2/3 + 7). 



PACS numbers: 12.15.Hh, 11.30.Er, 13.25.Hw 



4 



I. INTRODUCTION 

The Cabibbo-Kobayashi-Maskawa (CKM) quark 
flavor-mixing matrix [l[ provides an elegant explanation 
of the origin of CP violation within the Standard Model. 
CP violation manifests itself as a non-zero area of the 
unitarity triangle [2j. While it is sufficient to measure 
one of the angles to demonstrate the existence of CP vio- 
lation, the unitarity triangle needs to be over-constrained 
by experimental measurements in order to demonstrate 
that the CKM mechanism is the correct explanation of 
this phenomenon. Precision measurements of the sides 
and angles of the unitarity triangle are the focus of 
the physics program at the B Factories. While several 
theoretically clean measurements of the angle j3 exist Q , 
constraining the other two angles a and 7 is significantly 
more challenging. A theoretically clean measurement of 
sin(2/3 + 7) can be obtained from the study of the time 
evolution for B° -> D^~ir + and B° -> D^'p+ 
decays, which are available in large samples at the B 
factories, and for the corresponding CKM-suppressed 
modes B°->D^+tt- and B°->D^+p- [|. Measure- 
ments of CP asymmetries in decays B^-^D^^ir^ and 
S°^I? =F p ± decays have recently been published 

The interpretation of CP asymmetries in 
_B°^Z?(*) =F 7r ± decays as a measurement of sin(2/3 + 7) 
requires knowledge of the ratios of the decay amplitudes, 



r{D { *\) = 



A(B°^D(*)-%+) 



(1) 



However, direct measurements of the doubly Cabibbo 
suppressed branching fractions B(B°^D t -*' >+ Tr~) and 
B(B° — > D^ + p~) are not possible with the currently 
available data samples due to the presence of the copi- 
ous background from B°^L> ( *)+7r~, On the 
other hand, assuming SU(3) flavor symmetry, r(D^Tr) 
can be related to the branching fraction of the decay 
B ° _> Di* )+ ir- @: 



r (£>(*) vr) = tan(9 c 



/„(„ V B(B a -> DW-ttI 



(2) 



where 9 C is the Cabibbo angle, and fj^M / f d m is the ratio 

of flW and Dg meson decay constants 0, [^, flfj| . Other 
SU(3)-breaking effects are believed to affect r(D^n) by 
(10-15)% p. 

The dominant Feynman diagrams for the decays B° — > 



DW-7r+(p+), B° -> DW+7r-(p-), 5° 



(*)+ 



and £° -> £)i* } If M+ are shown in Fig. [TJ Since 
_B — > Ds k~ has four different quark flavors in the 
final state, a single amplitude contributes to the decay 
(Fig. [TJ;) . On the other hand, two diagrams contribute 
to B° -> L>M-7T+ and B° -> I>M+7r - : tree amplitudes 
(Fig. [JJi,b) and color-suppressed direct W^-exchange am- 
plitudes (Fig.[TJi,e). The latter are assumed to be negligi- 
bly small in Eq. ©. The decays B° -> D ( * ] ~ K+ (Fig.QJ) 
probe the size of the VF-exchange amplitudes relative to 
the dominant processes B° — > D^~tt + . 




^ (P ), 



FIG. 1: Dominant Feynman diagrams for (a) CKM-favored 
decays B° -> D ( * )_ 7r + (p + ), (b) doubly CKM-suppressed de- 
cays B° -> D ( * )+ 7r'(p"), and (c) the SU(3) flavor sym- 
metry related decays B° — > Z)i*' + 7r~ (p~); (d) the color- 
suppressed VK-exchange contributions to B° — » D'*'~7r + (p + ), 
(e) B° -> D ( * )+ 7r-(p-), and (f) decay B° -> A" w+ . 

The rate of 5° -> Di* )_ .ffW+ decays could be en- 
hanced by final state rescattering [k|, in addition to 
the VF-exchange amplitude. Such long-distance ef- 
fects could also affect the vector meson polarization 
in B° — > D*~ K* + decays. The angular distribution in 
vector- vector decays B° — > D* S V (V — p, K*) is given by 



'Deceased 

T Now at Temple University, Philadelphia, Pennsylvania 19122, 
USA 

*Now at Tel Aviv University, Tel Aviv, 69978, Israel 

§Also with Universita di Perugia, Dipartimento di Fisica, Perugia, 

Italy 

^Also with Universita di Sassari, Sassari, Italy 



d 2 r 



oc [(1- f L )(l + cos 2 9 D *J sin 2 6 V 



+ 4/ L sin 2 9 D , cos 2 By] , 



(3) 



where Qd* and dy are the helicity angles of D* + and the 
vector meson V, respectively, fa = |^4o| 2 /(^1^a| 2 ) is the 
longitudinal polarization fraction, and Aa=-i,o,+i are f ne 
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helicity amplitudes. These distributions are integrated 
over the angle between the decay planes of D* s + and V. 

For amplitudes dominated by the short-range (elec- 
troweak) currents, fi, is predicted to be near unity [l3l |. 
with corrections of order 0(my/m B ), where my is the 
mass of the vector meson produced by the weak current, 
and tub is the mass of the B meson. Thus, the measure- 
ment of Jl can constrain the size of the long-distance 
contributions in B° — ► D* s ~ K* + decays fl^ . 

The branching fractions B(B a — > £)« + 7T _ ) and 

B(B° — * Di*^ K + ) have been measured previously by 
the BABAR Collaboration In this paper we present 
the first evidence for the decays B° — > D* + p~ and 
B" -> Di*^K* + , and a limit on the rate of B° -> D+p". 
We also update the measurements of the branching frac- 
tions B(B a -► L>i* )+ 7r-) and B(B° ^ D { s *^ K+) with 
improved precision, using a 65% larger dataset. 



II. DATA SAMPLE AND THE DETECTOR 

We use a sample of 381 x 10 6 T(45) decays into BB 
pairs collected with the BABAR detector at the PEP-II 
asymmetric-energy e + e~ collider (l5| . A detailed de- 
scription of the BABAR detector is available elsewhere . 
The components of the detector crucial to this analysis 
are summarized below. 

Charged particle tracking is provided by a five-layer sil- 
icon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH). For charged-particle identification, ionization en- 
ergy loss (dE/dx) in the DCH and SVT, and Cherenkov 
radiation detected in a ring-imaging device (DIRC) are 
used. Photons and neutral pions are identified and 
measured using an electromagnetic calorimeter (EMC), 
which comprises 6580 thallium-doped Csl crystals. These 
systems are mounted inside a 1.5-Tesla solenoidal super- 
conducting magnet. We use the GEANT4 [17| software 
to simulate interactions of particles traversing the BABAR 
detector, taking into account the varying detector condi- 
tions and beam backgrounds. 



III. EVENT SELECTION AND ANALYSIS 

The selection of events of interest proceeds in two 
steps. First, we preselect events with at least three re- 
constructed charged-particle tracks and a total measured 
energy greater than 4.5 GeV, as determined using all 
charged particles and neutral particles with energy above 
30 MeV. In order to reject e + e~ — > qq(q = u, d, s, c) con- 
tinuum background, the ratio of the second to zeroth 
order Fox- Wolfram moments [l8j must be less than 0.5. 

Candidates for Df mesons are reconstructed in the 
L>+ -> 07T+, K° S K+ and K_*°K+ final states, with 
<j)->K + K-, K^tt+tt-, and K*°^K-tt+. The K° can- 
didates are reconstructed from two oppositely-charged 
tracks, and their momentum is required to make an angle 



I ^flight | < 11° with the line connecting their vertex and 
the e + e~ interaction point. All other tracks are required 
to originate from the e + e~ interaction region, loosely de- 
fined by | do | < 1.5 cm and \zo\ < 10 cm, where do and zo 
are the distances of closest approach to the primary e + e~ 
vertex in the directions perpendicular and parallel to the 
beams, respectively. In order to reject background from 
D + ^K^tt + or K*°tt + , the K + candidate in the recon- 
struction of Df^K^K + or K*°K + is required to satisfy 
positive kaon identification criteria, which have an effi- 
ciency of 85% and a 5% pion misidentification probabil- 
ity. The same selection is used to identify kaon daughters 
of the B° and K*+ mesons in decays B° -> D { * ] ~ . 
The selection is based on the ratios of likelihoods for 
kaon, pion, and proton identification in the SVT, DCH, 
and DIRC. The detector likelihoods are calibrated over 
a wide range of momenta using particles identified kinc- 
matically in clean decay chains, such as D* + — > D°ir + , 
D° — > K~ir + , and A — ► pir~ . In all other cases, kaons are 
not positively identified, but instead candidates passing 
a likelihood-based pion selection are rejected. The selec- 
tion efficiency of this "pion veto" is 95% for the kaons and 
20% for the pions. Pion daughters of B° and p~ mesons 
in the decays B° — > Ds + tt~ and B° — ► Di*^ + p~~ are re- 
quired to be positively identified. Decay products of 0, 
K*°, D+, and B° candidates are constrained to originate 
from a single vertex. 

We reconstruct p + — ► 7r + 7r° candidates by combining 
a well-identified charged pion with a ir° — > 77 candidate. 
The K* + candidates are reconstructed through the de- 
cays K*+ -> K+tt° and K*+ -> K° s n+ . The neutral 
pion candidates are reconstructed from a pair of pho- 
tons each with a minimum energy of 30 MeV. The in- 
variant mass of the photon pair is required to be within 
±25 MeV/c 2 of the nominal value 20]. The selected can- 
didates are constrained to the nominal ir° mass before 
forming the p + or K* + candidates. We require that the 
invariant mass of the two pions forming the p~ candidate 
be within ±320 MeV/c 2 of the nominal value [2(| , and the 
invariant mass of the K + ir° and KgTr + pairs be within 
±75 MeV/c 2 of the nominal K* + mass [loj]. ^5^+ pairs 
are constrained to a common geometric vertex. 

We reconstruct D* + candidates in the mode 
D* + ^Dfj by combining Df and photon candidates. 
Photon candidates are required to be consistent with 
an electromagnetic shower in the EMC, and to have an 
energy greater than 100 MeV in the laboratory frame. 
When forming a D* + , the candidate is required to 
have an invariant mass within 10 MeV/c 2 of the nominal 
value 20]. For B° -> D*+ p~ and B° ->■ D*~K* + modes, 
we apply a "7r° veto" by rejecting photons that in com- 
bination with any other photon in the event form an in- 
variant mass that falls within 125 < to 77 < 145 MeV/c 2 . 

The efficiency of the initial preselection dis- 
cussed above varies between 14% (B° — > D* + p~ , 
D+^K*°K+) and 48% {B° -»• D+%-, L>+ -> </>tt+). 
After the preselection, we identify signal B decay candi- 
dates using a likelihood ratio Rl = £ s i g /(£ s i g + £bkg)j 
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where £ S j g = Y[jV S i g (xk) is the multivariate likelihood 
for the signal hypothesis and £bkg = Yli Pbkgixk) is the 
likelihood for the background hypothesis. Here Xk repre- 
sents one of the discriminating variables described below, 
which are computed for each event. The likelihoods for 
the signal and background hypotheses are computed as 
a product of the probability density functions (PDFs) 
V s i g (xk) and Vhkg{xk), respectively, for the following 
selection variables: invariant masses of the cj>, K*°, p + , 
K* + , and K® candidates; x 2 confidence level of the 
vertex fit for the B° and Df mesons; the helicity angles 
of the (f), K*°, p + , K* + , and D* + meson decays; the 
mass difference Am(D* s + ) = m(D* s + ) - m(D+); the 
polar angle 9b of the B candidate momentum vector 
with respect to the beam axis in the e + e _ center-of-mass 
(cm.) frame; the angle 9t between the thrust axis of the 
B candidate and the thrust axis of all other particles in 
the event in the cm. frame; the event topology variable 
J 7 , and the kinematic variable AE, described below. 
Correlations among these variables are small. 

The helicity angle 9h is defined as the angle between 
one of the decay products of a vector meson and the 
flight direction of its parent particle in the meson's rest 
frame. Polarization of the vector mesons in the signal de- 
cays causes the cosines of their helicity angles to be dis- 
tributed as cos 2 8 H (0, K*°, p + , and K*+) or 1 - cos 2 9 H 
(D* + ), while the random background combinations tend 
to produce a more uniform distribution in cos 9h, with 
a peak in the forward direction (which corresponds to 
a low-energy ir°) for p + and K* + candidates. We do 
not include the helicity angles for D* + , p + , and K* + 
mesons in the likelihood ratio Rl for the vector-vector 
B" -> D*+p~ and B° -> D*~ K*+ modes, since the po- 
larizations of the vector mesons in these decays are not 
known a priori. Instead, the helicity angles are used in 
the multi-dimensional likelihood fit to determine the po- 
larizations, as discussed below. 

The variables cos 6b, cos 9t, and T discriminate be- 
tween spherically-symmetric BB events and jet-like con- 
tinuum background using event topology. The polar an- 
gle 9b is distributed as sin 2 9b for real B decays, while 
being nearly flat in cos 9b for the continuum. BB pairs 
form a nearly uniform |cos#t| distribution, while the 
| cos 9t I distribution for continuum events peaks at 1. A 
linear (Fisher) discriminant T is derived from the values 
of sphericity and thrust for the event, and the two Leg- 
endre moments Lq and L2 of the energy flow around the 
.B-candidate thrust axis [l9l |. 

The ratio Rl has a maximum at Rl — 1 for sig- 
nal events, and at Rl = for background originat- 
ing from continuum events. It also discriminates well 
against B decays without a real -D+ meson in the final 
state. The Monte Carlo (MC) simulated distributions 
of the Rl variable for signal and background events, in 
B° — > Dfir~ decays, are shown in Fig. O 

Finally, two other variables Toes and AE take advan- 
tage of the unique kinematic properties of the e + e~ — > 
T(4S) — > BB decays. The beam energy spread is signif- 




FIG. 2: Distribution of the likelihood ratio Rl for the mode 
B° -> D+tt~, D+ -> 07T+. Shown are (a) the simulated signal 
events, and (b) the sum of the simulated background samples 
from the B° and B + decays, and e + e~ — > qq events. 



icantly smaller than the energy resolution of the recon- 
structed B mesons, and at the same time larger than the 
momentum resolution. The momentum of the signal can- 
didates is included in the beam-energy-substituted mass 
toes = y/(s/2 + Pi -PB) 2 /E'f -p B , where yfs is the to- 
tal cm. energy, (£^,pj) is the four-momentum of the 
initial e + e~ system, and p_e is the B° candidate momen- 
tum, both measured in the laboratory frame. The second 
variable is AE = E* B - y/s/2, where E* B is the B° candi- 
date energy in the cm. frame. For signal events, the Toes 
distribution is nearly Gaussian and centered at the B 
meson mass with a resolution of about (2.5-2.8) MeV/c 2 , 
and the AE distribution has a maximum near zero with a 
resolution of (17-25) MeV. We include AE in the defini- 
tion of the likelihood ratio Rl ; toes is used as a discrim- 
inating variable in the maximum likelihood fit described 
below. 

We parameterize the signal and background PDFs 
using large samples of simulated events. We select 
B° -> d{* )+ tt- and B° -> K+ candidates that 

satisfy Rl > 0.85, and accept B° — ► p~ and 

B° -> L>i* )_ /f*+ candidates with R L > 0.96. We mea- 
sure the relative efficiencies er l of the Rl selections in 
copious data samples of decays B° — > D~u + . D~p + 
(D~ -> K+7C-TT-, Kjn-) and B+ -> D*°tt+, D*°p+ 
(D*° -> D°7, D° -> K-JT+) in which the kine- 
matics is similar to those of our signal events, and 
find that they are consistent with Monte Carlo esti- 
mates of e Rl ps 70% (40%) for B° -> D { * )+ Tt- and 
B° D ( ; ] -K+ (B° -> D ( ; )+ p- and B° -» D ( s * y K*+) 
modes. The fraction of continuum background events 
passing the selection varies between 2% and 15%, de- 
pending on the mode. 

Less than 30% of the selected events in the 
B°^D S *V, B° ->■ D*~K+, B Q ^Di* )+ p-, and 

B° -> D { * ] ~K*+ channels (< 2% in B° -> D+tt- and 
B° — > D~K + ) contain two or more candidates that sat- 
isfy the criteria listed above. In such events we select a 
single B candidate based on an event \ 2 formed from 
AE, m(D s ) and (where appropriate) Am(D* + ), m p , 
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rriK*, m^o and rnxs, and their uncertainties. 



IV. EXTRACTION OF SIGNAL YIELDS 

After the Rl requirement is applied, we define the 
region of interest using the beam-energy-substituted 
mass ibes and the mass of the Dj candidate m(D s ). 
We require 5.2 < to E s < 5.3 GeV/c 2 and \m(D s ) - 
m(D s ) PDG \ < 50 MeV/c 2 for B° -> Dfir-, B° -> D+p~, 
and B° — > D~A^*) + modes, where m(-D s )pDG is the 
world average D s mass [20]. The invariant mass m(D s ) 
has a resolution of (5-6) MeV/c 2 , depending on the 
decay mode. The selection is significantly broader than 
the region populated by the signal events, and allows 
us to constrain backgrounds in the signal region. For 
B° -> D* s +tt-, B° -► D*+p~, and B° -► we 
require |m(D s ) - m(D s ) PDG | < 10 MeV/c 2 . 

Five classes of background events contribute to the 
fit region. First is the combinatorial background, in 
which a true or fake Z?s candidate is combined with 
a randomly-selected light meson. Second, B meson 
decays such as_B°^D ( * )+ 7r" or B°^D^+p- with 
D + ^K®,ir + or A*°7r + can constitute a background for 

the B° -> Di* )+ 7r" and £° -> £>£* )+ p- modes if the pion 
in the D decay is misidentified as a kaon (reflection back- 
ground). The reflection background has nearly the same 
itjes distribution as the signal but different distributions 
in AE and m(D 8 ). The corresponding backgrounds for 
the B° -> D~K(*>+ mode (B°^D- K^+) are negligi- 
ble. Third, rare charmless B decays into the same fi- 
nal state, such as B°—>K { * )0 K + h (where h = it, p, K, 
or K*), have the same toes and Ai5 distributions as 
the B° — > _D s /i signal, but are nearly flat in m(D s ). 
The charmless background is significant in B — * D+ir~ , 
B° -> D+p-, and £° -► D~ decays, but is effec- 
tively rejected by the Am(D* + ) variable for the modes 
with D*+. 

Finally, two classes of background events have nearly 
the same distribution as the signal events in both m(D s ) 



and toes- Fot B° — > Z?i*^ K* + modes we take into ac- 
count the potential contributions from the non-resonant 
decays B° — > D s K°ir + (which have recently been ob- 
served and color-suppressed B° — > fi*' ~A" + 7r° 
(unobserved so far). Analogous non- resonant modes 
5° — * Z^i*^^ - ^ require the additional popping of 
a color-matched pair. They are expected to be 
small compared to B a — > D^^p - [H[ and are ignored. 
Finally, crossfeed background from misidentification of 
_B° -> £>£* )- 7r+ events as 5° -> D { * ] ~ K+ signal, and vice 
versa, needs to be taken into account. 

For each mode of interest, we perform an unbinned 
extended maximum-likelihood (ML) fit to separate the 
signal events from the backgrounds and extract the sig- 
nal branching fractions. For B° — > Dfir~, B° — > Dfp~, 
B° — > D~ K + , and £>° — > D~ K* + , we perform a two- 
dimensional fit to the Toes and m(D s ) distributions. For 



B° -> D*+7T" and 5° -> D*- K+ decays, we fit the one- 
dimensional Toes distribution. In vector-vector modes 
B° -> L>*+p"and 5° -> D*~ A*+, we constrain both the 
branching fractions of the signal modes and the po- 
larization of the vector mesons by performing a three- 
dimensional fit to the distribution of toes, and the two 
helicity angles of the D* + and p~ (K* + ) mesons. 

For each B decay, we simultaneously fit distributions 
in the three Df decay modes, constraining the signal 
branching fractions to a common value. The likelihood 
function contains the contributions of the signal and the 
five background components discussed above. The func- 
tion to be maximized is 



exp 



n kr , 



n 

i=l 



(4) 



where rij m is the number of events for each event type j 
(signal and all background modes) in each D s decay mode 
to, and VjmiCi) is the probability density function of the 
variables Q = (toes, m(D s ), cos 9d% , cos 9y) for the zth 
event. The likelihood product is computed over all can- 
didates A can d in the region of interest. We parameterize 
the event yields as 



N BS BiB. 



where to stands for Di 



D+ 



K°K+ 



N 



BB 



4>TT + 

381 x 10 6 



'm 7 (5) 

D+ -> K*°K+, or 
Bj is the B decay 



branching fraction, B^ s is the branching fraction of the 
TO-th mode, and e m is the reconstruction efficiency. 

The branching fractions of the channels contributing 
to the reflection background are fixed in the fit to the 
current world average values [2(| and the branching frac- 
tions of the crossfeed backgrounds are determined by it- 
erating the fits over each B decay mode. The branch- 
ing fractions of the non-resonant backgrounds are fixed 
to the values recently measured by BABAR [2l|. In 



the case of B° -> D 
to B° -» D { ; ] -K*+ ( 
branching fraction by 



(*)- 



A + 7T°, which can contribute 
— » A + 7t°), we estimate the 



B(B° -» D^-K+Ti a ) m 
B{B+ -> D^-K + n+) 



B(B°^D°7t") 
B(B+->D°tt+) 



(6) 



This scaling assumes that the dominant mechanism for 
producing both D^*' K + tt° and D{*^ K + ir + final states 
is a sub-threshold production of a charmed D**° meson, 
which subsequently decays into D s *^ K + , as indicated 
by the invariant mass spectrum of D s K + [2lj . Since 
the decay B° — > D**°tt° is color-suppressed compared to 
B + — > D**°ir + , we estimate the color suppression factor 
from the B° — > D°ir Q decays. Direct production of the 
color-suppressed Di*^ K + tt° final state (without the in- 
termediate D**°) results in a smaller branching fraction 
estimate. We assign a 100% systematic uncertainty to 



B(B° 



D 



(*)- 



K+ir° 
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The expected yields of the dominant B-decay back- 
grounds are listed in Table HI 

The PDFs and efficiencies for the signal, reflec- 
tion, and crossfeed backgrounds are determined inde- 
pendently for each decay mode using Monte Carlo 
samples. The signal contribution is modeled as a 
Gaussian (B° -> Dfn~ and B Q -> D~K + ) or a "Crys- 
tal Ball" function [22] in ttjes and a double Gaus- 
sian in m(D s ). The combinatorial background is de- 
scribed in rriES by a threshold function [231 ] . dN/dx oc 
xy/l — 2x 2 /sexp [— £ (l — 2x 2 /s)] , characterized by the 
shape parameter £. This shape parameter, common to 
all Df modes, is allowed to vary in the fit. In m(D s ), the 
combinatorial background is well described by a combi- 
nation of a first-order polynomial (fake £>+ candidates) 
and a Gaussian with (5-6) MeV/c 2 resolution (true Df 
candidates). The charmless background is parameterized 
by the signal Gaussian shape in thes & n d a first order 
polynomial in m(D s ). 

Ideally, the distribution of the helicity angles in the 
vector- vector decays is given by Eq. ([3]) . The helicity an- 
gle 9d* is defined as the angle between the direction of 
the photon in D* — > D s j and the direction of the B in the 
rest frame of the D* s candidate. The helicity angle 9y is 
similarly defined by the direction of the charged daughter 
particle in the decays p + — ► tt + tt°, K* + — ► K + ir°, and 
K* + — ► A^7r + . Since the momenta of the decay products 
in the laboratory frame depend on the helicity angles, ac- 
ceptance and efficiency effects modify the ideal angular 
distribution. We determine the PDFs of the signal events 
using the Monte Carlo simulation, and measure the an- 
gular distribution of the combinatorial background in the 
data region rriES < 5.27 GeV/c 2 . 

For B° -> D+TT-, B" -> D+p-, and B° -> DjK+, the 
fit constrains 14 free parameters: the shape parameter 
of the combinatorial background £ (1 parameter for all 
D+ modes), the slopes of the combinatorial and charm- 
less backgrounds in m(D s ) (3 parameters), the fractions 
of true candidates in combinatorial background (3), 
the numbers of combinatorial background events (3), the 
numbers of charmless events (3), and the branching frac- 
tion of the signal mode (I). In the B° — > D~ K* + mode 
(6 individual sub-modes, spanning 3 channels and 2 
K* + channels), 26 free parameters are constrained. For 
the B° -> D* s +tt- and B° -> D*~ K+ decays, 5 free pa- 
rameters are determined by the fit: £ (1 parameter for 
all Df modes) , the number of combinatorial background 
events (3), and the branching fraction of the signal mode 
(1). For 5° -> D*+p- and B° -> D*-R*+ fits, we add 
one more free parameter to the fit: the longitudinal polar- 
ization fraction (see Eq. The total number of free 
parameters is 6 in B° -> D*+p- and 9 in B° -> D*~K* + . 

The results of the fits are shown in Figs. OH] and sum- 
marized in Table HT1 



V. SYSTEMATIC UNCERTAINTIES 

For the branching fractions, the systematic errors are 
dominated by the 13% relative uncertainty for B(Df — > 
(f>n + ) 20]. The uncertainty in the branching fraction 
ratio B(D+^K*°K+)/B(D+^(t>n+) contributes (2-4)%, 
depending on the decay channel. For B(D£—>K® K + ), we 
use the most recent measurement from the CLEO Collab- 
oration (24|, which differs from the previously reported 
central value [2(| by about 50%. We estimate uncertain- 
ties due to modeling of the resonance (K*°, <j>, p, and 
K* + ) lincshapes by measuring the effect of the lineshape 
variation on signal selection efficiency. 

The uncertainties in the signal selection efficiency are 
determined by the accuracy with which the detector ef- 
fects are modeled in the Monte Carlo simulations. Track- 
ing, particle identification (PID), photon, ir° and Kg 
reconstruction efficiencies are measured across the wide 
range of particle momenta in the dedicated data control 
samples. The tracking efficiency and resolution are ad- 
equately reproduced by the simulations. The simulated 
distributions are corrected for the efficiency and resolu- 
tion of the 7T° reconstruction. The efficiency of the Rl 
cut is also measured in the data control samples, as dis- 
cussed in Section IIIII 

The uncertainties due to the knowledge of the signal 
and background PDFs in the ML fit are estimated by 
measuring the variation of the fitted values of the branch- 
ing fractions when PDF parameters are varied within 
their uncertainties. The correlations between parame- 
ters are taken into account. The uncertainties in the sig- 
nal PDF parameters for the key discriminants AE, jbes, 
m(D s ), Am(Dl + ), and cos9 D *+ are determined by com- 
paring data and Monte Carlo simulations for the samples 
of decays D~tt+, D~ p+ (£>~ K+tt-tt-, K°Tr~) 

and B+ -> D*°n+, D*° p+ (D*° -> L>° 7 , D° -> K~-K + ). 
The uncertainties in the signal PDFs for cos# Pj i<> and 
the PDFs for the peaking backgrounds arc determined 
by Monte Carlo simulations. These distributions depend 
on the modeling of the charged track and ir° reconstruc- 
tion, discussed above. The helicity angle PDFs for the 
continuum background are determined in the data side- 
band ities < 5.27 GeV/c 2 , and their uncertainties are 
statistical in nature. 

Uncertainties due to reflection and crossfeed back- 
grounds include the uncertainties in the branching frac- 
tions of the relevant modes, and also account for the 
contributions of the sub-dominant modes that are not 
explicitly included in the ML fit. These contributions 
dominate the systematic uncertainty for the B —> p~ 
mode, which has a small absolute branching fraction. 

As ML estimators may be biased in small samples, 
we measure the bias using a large ensemble of simu- 
lated experiments. In each of these experiments, gen- 
erated according to the sample composition observed 
in data, the signal and _B-decay background events 
are fully simulated, and the combinatorial background 
events are generated from their PDFs. The bias is 
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found to be negligible for the 1- and 2-dimensional ML 
fits {B° Di* )+ ir-, B°^D ( s t) -K+, B°^D-K*+ 
modes). On the other hand, we find that in the vector- 
vector modes (B° -> D*+ p~ and B Q -> D*~ K*+ de- 
cays), the 3-dimensional ML fits underestimate the true 
values of the signal branching fraction and the fraction of 
the longitudinal polarization. We correct for the biases 
of Ai3 = (-0.37±0.03) x 10" 5 and Af L = (-5.3±0.6)% 
(B° -> D*+p-) and AB = (-0.14 ± 0.04) x 10~ 5 and 
A/l = (-5.5 ± 0.8)% (B° -> D*~K*+). We assign a 
conservative uncertainty of 50% of the bias to this cor- 
rection. 

For the longitudinal polarization fractions in the 
vector-vector modes, the systematic errors are dominated 
by the uncertainties in the shapes of the signal and back- 
ground PDFs and the fit bias. The systematic uncertain- 
ties for each mode are summarized in Tables IHUfVl 



VI. RESULTS 

We estimate the significance of a non-zero signal yield 
by computing S = y/-2 log(£ /£max), where £ max is 
the maximum likelihood value, and Co is the likelihood 
for a fit in which the signal contribution is set to zero. In- 
cluding systematic uncertainties and assuming Gaussian- 
distributed errors, we obtain signal observation signifi- 
cances of 3.9 (B° -> D* s +p-), 4.6 (B° -> D~K* + ), and 
3.1 (B° — > D*~K* + ) standard deviations, providing the 
first evidence for these decays. We test that S measures 
the probability for the background events to fluctuate to 
the observed number of signal events with a large ensem- 
ble of simulated experiments. For each such experiment, 
we generate a set of pure background events according 
to the PDFs and sample composition observed in our 
dataset. We then fit each simulated experiment and mea- 
sure the signal and background yields and, for the vector- 
vector modes, the polarization fraction f^. By counting 
the fraction of such pseudo-experiments in which the sig- 
nal yields are at least as large as the yield observed in 
the real dataset, we confirm that S 2 follows closely the 
X 2 distribution with one degree of freedom. 

The branching fraction and polarization results are col- 
lected in Table [TTJ Since we do not observe a significant 
event yield in B° — ► Dfp~, we set a 90% confidence- 
level Bayesian upper limit assuming a constant prior for 
B(B° D+p-) > 0. 



states B° -> D*+p~ and B° -> D*-R*+: 
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= [4.4_ 1 2 ± U.oJ X 1U 


£ t 7">0 

/m-B - 


D *s + P~) 


= 0.8oT 28 ± 0.15 


B(B U - 


-> D~K+) 


To n i in a \ d ol -. , ~i n — 5 

= [2.9 ± 0.4 ± 0.3J x 10 


B(B° - 


> D*~K+) 


= [2.4 ±0.4 ±0.2] x 10~ 5 


B(B° - 


> DjK*+) 


= [3.6±J;g ± 0.4] x 10~ 5 


B(B° -> 


D*-K*+) 


= [3.0±\i ± 0.3] x 10~ 5 


h(B° -> 


D*-K*+) 


= Q.96±g:|f ± 0.08 , 



where the first quoted uncertainty is statistical, and the 
second is systematic. 

The branching fractions for B° — > Di*'~ arc 
small compared to the dominant decays B° — > Z)(*)~7r + 
and S° — > D^*' p + , implying insignificant contri- 
butions from the color-suppressed VF-exchange dia- 
grams. The ratios £(5° -> D-K+)/B{B° -> D*- K+) 
and 6(B° -> D-K*+)/B{B° -> D*~K*^) are consis- 
tent with the expectation of unity [12|. The pre- 
dictions for the branching fractions of B u ^ D s k~ 
and -B — * Z?s + p - decays are based on the fac- 
torization assumption [25j and depend on the es- 
timates of the hadronic form factors. The ob- 
served pattern B(B° -► .D+p") < S(B° D+tt~) w 
6(5° D*^tt~ ) < 6(5° D*+p-) appears to be most 
consistent with the form factors computed in 26]. The 
polarizations of the vector mesons in B° — > D* + p~ and 
B° — > D*~ K* + are consistent with expectations [T3. fl3j . 

Assuming the SU(3) relation, Eq. and the recent 
value of j '„(») / j £,(,) from an unquenched Lattice QCD 
calculation Q, we determine the ratios of the CKM- 
suppressed to CKM-favored decay amplitudes in decays 
B° -> DW ± tt t and B° D^p^: 

r(Dir) = [1.75 ± 0.14 (stat) ± 0.09 (syst) ± 0.10 (th)]% 

r{D*7r) = [1.81±g;^ (stat) ± 0.12 (syst) ± 0.10 (th)]% 

r{Dp) = [0.7lt^26 (stat) ±0.11 (syst) ± 0.04 (th)]% 

r{D*p) = [1.50±g;|? (stat) ± 0.16 (syst) ± 0.08 (th)]% 



VII. CONCLUSIONS 

We report the following improved measurements of the 
branching fractions for the rare decays B° — * 0^ + 7r" 
and B° — * Z?i*^ K + , and the first measurements of the 
branching fractions for the decays B° — > p~ and 

B a — > Ds*'' as well as the measurements of the lon- 

gitudinal polarization fractions Jl in vector-vector final 



where the first error is statistical, the second includes 
experimental systematics, and the last accounts for the 
uncertainty in the theoretical value of /r>(*)//D(») @- 
These amplitude ratios are below 2%, which implies small 
CP asymmetries in B°^D^w ± and B°^D^p ± 
decays, making it difficult to measure sin(2/3 + 7) 
precisely in these decays. The results presented here 
supersede our previously published measurements [lij]. 
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TABLE I: Expected background yields from the dominant B decay modes, fixed in the likelihood fit. 



Sig 


nal mode 






Background mode 




N(D+ 


) N(Df -> K* a 


K+) N(D+ -> K° S K+) 


1st 


* U s 7T 








D7if+ 




1 c i n o 
1.0 It U.Z 


0/1-1-0 1 
U.^ zt U. 1 


n q 4- o o 

U.O it U.U 
















171 +1 9 


91 1 + 1 fi 


i q q _|_ i n 
±o.y it ±.u 


B 









B° 


-£>:-*-+ 




0.9 ± 0.2 


0.3 ± 0.1 


0.3 ± 0.0 












-> D~p+ 




0.5 ± 0.1 


3.5 ± 0.7 


1.8 ± 0.4 










7->0 




— ► i7 71 


U.O m U. 1 


l.Z It u.z 


O fi 4- O 1 
U.o it U.l 


B 








B° 






6.9 ± 2.0 


1.4 ± 0.4 


1.6 ± 0.4 










B° 






fi ^ -I- 1 1 
0.0 i 1.1 


i q _|_ n 9 
1.0 It XJ.Z 


J..D it U.O 












r\*4- 




o.y it u. j 


11-1-09 

± . _L HI U.Zr 


i o 4- n 9 

_L .U It U.Z 










JD 


^ 71 ~ 

— > .L/ p 




26.0 ± 4.4 


35 2 ± 5 9 


30 1 ± 5 










JD 


— > i> p ,L> 




0.3 ± 0.0 


6.1 ± 3.7 


8.5 ± 1.4 










B° 


— > p , u 


— ► D~7T° 


0.9 ±0.4 


1.5 ±0.6 


2.2 ±0.5 


B° 








rO 


— > JJ p 




0.7 ± 0.2 


1.7 ± 0.4 


2.6 ± 1.2 










R° 


— > _L» p ,.D 


' -» D-ir° 


0.1 ± 0.0 


0.8 ± 0.1 


0.8 ±0.1 


11) 










Dtn- 




U.O IE U. 1 


U.O it U.U 


9 4-0 
U.Z It U.U 












->D*-K+ 




i 4 _|_ n q 

l.^t It U.O 


n q _|_ n 1 

U.O It U. 1 


n 9 4- o o 
u.z in u.u 










7->0 


— * D~K + 




n n _i_ n q 


9 9 4-no 
z.z it \j.z< 


i q _i_ n 4 
1.0 It U. 1 ! 


B 


— 

— > D* K + 






B° 






0.9 ± 0.1 


0.2 ± 0.0 


0.1 ± 0.0 










S° 






10-1-00 

i.u it u.z 


0/14-01 
U.4 it U.l 


u.o it u.i 










B° 


o;+7T- 




n c _(_ n i 
U.O it U. 1 


9 4-00 
U.Z it U.U 


9 4-0 
U.Z it U.U 


B" 


-> D~K*+ 


, K* + 


-I- 


B° 






0.4 ± 0.2 


0.1 ± 0.0 


0.1 ± 0.0 










B° 






1.9 ± 0.7 


0.8 ± 0.1 


0.6 ±0.1 


B° 


D~K*+ 


, K*+ 




B° 






2.6 ± 0.4 


0.9 ±0.2 


0.9 ±0.2 










B° 






1.1 ±0.5 


0.3 ±0.2 


0.4 ±0.2 










B° 






0.4 ±0.4 


0.1 ±0.1 


0.1 ±0.1 


B° 


-> D*~K* 


f , 




B° 


-> L^ir^ 




0.2 ±0.0 


0.1 ±0.0 


0.0 ±0.0 










B° 


-> D*--k+K° 




0.6 ±0.4 


0.3 ±0.2 


0.2 ±0.1 


B° 


D*~K* 




h - if +7r° 


B° 






0.5 ±0.1 


0.1 ±0.0 


0.2 ± 0.1 










B° 


-> D*-K+tt" 




0.1 ±0.1 


0.0 ±0.0 


0.0 ±0.0 
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TABLE II: The number of reconstructed candidates (iV raw ), the signal yield (N B i g ), computed from the fitted branching fractions, 
the combinatorial background (N coul b), and the sum of charmless, reflection, non-resonant, and crossfeed contributions (-/V pea k), 
extracted from the likelihood fit. Also given are the reconstruction efficiency (e), the signal significance S, the measured 
branching fraction £>, and the fraction of longitudinal polarization /l (where appropriate). The first uncertainty is statistical, 
and the second is systematic. 



B mode 


D s mode 


raw 


iVsig 


comb 




s(%) 


S 


6(10~ b ) 


h 








582 


51 ±10 


500 ± 24 


32 ± 10 


25.2 








5° 


— > D + 7T~ 


Df^K*»K+ 


402 


19 ±4 


352 ± 20 


36 ±8 


8.0 


8.2a 


2.5 ±0.4 ±0.2 










282 


19 ±4 


245 ± 16 


25 ±7 


19.4 
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34 ±6 
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1.7 ±0.3 


16.7 








B° 
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r s /I 


D+^K*°K+ 


96 
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D+^K° S K+ 


52 
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41 ±7 
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B° 


— » F) + n~ 


D+^K*°K+ 
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1.3a 
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D+^K° 8 K+ 
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4±4 


544 ± 25 


70 ±13 


8.3 




< 2.4 (90% C.L.) 
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22 ±6 
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0.8 ±0.2 


6.3 








B° 


— > D*+n _ 

^ u s P 


D+^K*°K+ 
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1.9 


3.9cr 


4A±\ l ± 0.8 


°- 86 -0.28 ±0.15 






D+^K° S K+ 
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307 
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15 ±7 
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76 
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47 ±8 
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Df^K° s K+ 


34 
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21 ±5 
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95 
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K* 






45 


4± 1 


40 ±7 


1 ±2 
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Df^K° s K+ 


33 


3±1 


27 ±6 


1±3 
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9±3 
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1±4 
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K* 
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3±1 
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6 ±4 
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3±1 
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B u 
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0.8 ±0.4 
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K* 




D+^K*"K+ 
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2 ± 1 
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0.4 ±0.2 


2.5 












Df^K° s K+ 
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1 ± 1 


6±3 


0.2 ±0.1 
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B° 
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22 ±6 


0.6 ±0.2 
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3.1a 
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11 
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TABLE III: Relative systematic uncertainties for the branching fractions of B° — > D S X modes (%). 







B" -> D-R+ 


B U - Dfp- 


B" -> D~K* + 


^BB 


1.1 


1.1 


1.1 


1.1 


Tracking efficiency 


1.6 


1.3 


1.7 


2.0 


PID efficiency 


1.0 


1.0 


1.0 


1.0 


7T° efficiency 






3.0 


3.0 


Rl cut efficiency 


1.6 


1.4 


2.8 


2.0 


MC statistics 


0.8 


0.7 


1.8 


1.1 


Kg efficiency 


0.4 


0.3 


0.9 


0.0 


PDF parameters 


0.8 


0.7 


2.8 


1.4 


AE, m(D s ) PDFs 


0.6 


1.7 


3.2 


1.7 


B(Df -> 07T+) 


7.7 


8.3 


14.7 


7.7 


B(Df -> K*°K + ) 
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3.1 


1.8 


3.1 


B(D+ -> K° S K+) 


1.6 


1.4 


3.7 


0.6 


Reflection background 


2.0 


0.7 


10.2 


0.6 


Crossfeed background 


0.9 


0.4 


15.6 


2.5 


Resonant lineshape 


1.6 


1.7 


1.8 


4.9 


TOTAL 


9.3 


9.7 


25.0 


13.1 
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TABLE IV: Relative systematic uncertainties for the branching fractions of B° — > D* S X modes (%). 







B l > ^ D*-K + B u 


- D* 8 + p- 


B" -> D*-K*+ 


N BB 


1.1 


1.1 


1.1 


1.1 


Tracking efficiency 


1.6 


1.3 


1.7 


2.0 


PID efficiency 


1.0 


1.0 


1.0 


1.0 


7T° efficiency 


- 


- 


3.0 


3.0 


7 efficiency 


1.8 


1.8 


1.8 


1.8 


B{D*+ -> D s + 7 ) 


0.7 


0.7 


0.7 


0.7 


Am(D s « + ) 


1.0 


1.0 


1.0 


1.0 


cos D *+ 


1.0 


1.0 


1.0 


1.0 


tv° veto efficiency 


- 


- 


2.0 


2.0 


Fit bias 


- 


- 


4.1 


2.3 


Rl cut efficiency 


1.5 


1.7 


2.3 


1.6 
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0.8 


0.8 


1.1 


1.6 


Kg efficiency 


0.4 


0.0 
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0.3 
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1.9 


2.1 


3.0 


3.1 
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3.0 
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8.4 
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5.9 


B(L>+ 


2.7 


3.8 


3.4 


4.1 




1.5 


1.3 


0.9 


1.3 


Reflection background 


3.4 


0.0 


2.1 


0.0 


Crossfced background 


1.9 


0.4 


0.0 


1.6 


Resonant lineshape 


0.8 


0.4 


1.6 


0.3 


TOTAL 


10.5 


10.2 


18.2 


10.6 



TABLE V: Absolute systematic uncertainties for the longitudinal polarization fraction (%). 





B u - D'+p- 


B y) -> D*~K*+ 


Fit bias 


2.7 
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Rl cut efficiency 
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MC statistics 


0.6 


0.7 


K% efficiency 


0.4 


0.0 
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0.0 


0.5 


B(Df 4>n + ) 


0.8 


0.8 


B(D+ -> K*°K+) 
B(D+ -» K° S K+) 


0.7 


0.2 


0.5 


0.2 


Reflection background 


0.5 


0.0 


Crossfeed background 


0.0 


0.6 


Resonant lineshape 


0.6 


0.1 


TOTAL 


15.0 
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FIG. 3: (a,c,e,g) rriES projection of the fit with \m(Df) — m(D^")pDG| < 10 MeV/c 2 and (b,d,f,h) m(D s ) projection with 
5.275 < m ES < 5.285 GeV/c 2 for (a,b) B° -> D+tt", (c,d) B° -> DJ K+ , (e,f) 5° D+p~, and (g,h) B° -> DJ K*+ . The 
black solid curves correspond to the full PDF from the combined fit to all decay modes. Individual contributions are shown 
as solid red (signal), green dashed (combinatorial background), and blue dotted (sum of reflection, charmless, crossfeed, and 
non-resonant backgrounds) curves. 
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t — 1 — 1 — 1 — r 




FIG. 4: mES projection of the fit for (a) B° — » _Ds + 7r~and (b) B° — > D* s ~ K + . The black solid curves correspond to the full 
PDF from the combined fit to all decay modes. Individual contributions are shown as solid red (signal), green dashed 
(combinatorial background), and blue dotted (sum of reflection, charmless, and crossfeed backgrounds) curves. 
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FIG. 5: Projections of the fit on to (a,b) mES, ( c ) cos# p , (d) cosOk*, and (e,f) cosOd* for (a,c,e) B° — > Dg + p and (b,d,f) 
B a — > D*~ K* + . For helicity projections, a selection 5.275 < ules < 5.285 GeV/c 2 is applied. The black solid curves correspond 
to the full PDF from the combined fit to all decay modes. Individual contributions are shown as solid red (signal), green 
dashed (combinatorial background), and blue dotted (sum of reflection, charmless, crossfeed, and non-resonant backgrounds) 



